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Abstract 

Two Magma functions are given: one computes linear systems of plane 
curves with non-ordinary singularities and the other computes a scheme 
which parametrizes given degree plane curves with given singularities. 
These functions provide an efficient tool to construct explicit equations of 
singular plane algebraic curves. 

By computing singular branch curves, we obtain equations of normal 
quartic surfaces in CP"^ having the following combinations of rational dou- 
ble points: D5E7E7, DrDeDe, EeDgDs, EgDis, EeEeDr, EeEgDs, EyDeDe, 
E7D12, EyEeDg. These are all possible cases with total Milnor number 19 
which have no point of type A„. 
2000 MSG: 14Q05, 14H50. 



1 Introduction 

A quartic surface in CP'^ with only rational double points as singularities is 
an example of a K3 surface (Calabi-Yau variety of complex dimension 2), an 
important class of surfaces with applications to theoretical physics, in particular 
to string theory. Mathematicians have been studying K3 surfaces for over one 
hundred years, at least since [Hud]. 

Several authors have studied normal quartic surfaces with rational double 
points. A Dynkin graph of such a double point is of type A„, or Ep, p — 6,7 
or 8, with Milnor number the index n, m or p. Let r be the total Milnor number, 
i.e. the sum of all Milnor numbers. The maximum value for r is 19. Based on 
Nikulin's ([M]) and Urabe's f fUrl] , [Ur2| ) work, Yang (fYa^) has computed, for 
r = 19, 18 or 17, all possible configurations of rational double points that occur 
on a normal projective quartic surface in CP"^. 

The number of possibilities for r = 19 is 278, but, to my knowledge, equa- 
tions of such singular surfaces have not been given (Urabe's method is not 
constructive), except for the case of a quartic surface with an A19 singularity 



constructed by Kato and Naruki ( |KN] ) . Surprisingly, removing from Yang's list 
the cases which contain at least one singularity of type A„, only the following 
9 cases remain: D5E7E7, DrDgDe, EgDsDs, EgDis, EgEgDy, EgEsDs, ErDgDe, 
E7D12, EyEgDg. In this paper we compute an equation for each of these cases. 
We notice that our method can be used to construct other surfaces in the list. 
For instance an equation for the case A17A1A1 is not difficult to obtain. 

Let X be a quartic surface with a double point p. The projection from p 
gives X as a double covering of the plane. Such a covering is determined by its 
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branch locus B (the projection of the ramification curve to P^). After computing 
the defining equation of it is not difficult to obtain the equation of X. 

In this paper we give an efficient tool to compute singular plane curves. 
All computations are implemented with the Computational Algebra System 
MAGMA f |BCP| ). The Magma function LinearSystem computes only linear 
systems of curves with ordinary singularities. To overcome this restriction, we 
define the function LinSys which calculates systems of curves having (any type 
of) non-ordinary singularities. 

Let L be a linear system of plane curves of degree d. To obtain a curve of L 
with given singularities, one imposes conditions to its elements. If the number of 
conditions is greater than the dimension of L, this curve may not exist. Suppose 
that the function LinSys returns no sections (or returns a non-reduced curve). 
In this case the singular points pi, . . . ,pn (possibly infinitely near) may be in a 
special position. We give a function ParSch whose output is a scheme which 
parametrizes plane curves (or linear systems of curves) with given singularities 
at pi, . . .,Pn- 

These two functions can be used to construct any type of singular plane 
curves, depending the success mostly on computer power. 

The paper is organized as follows. In section [5] we introduce the functions 
LinSys and ParSch. The corresponding code lines, implemented with the Com- 
putational Algebra System Magma ([BCP]), are given in the Appendix. In Sec- 
tion [3] we relate the quartic surface X with its branch locus B and we present a 
table with possibilities for the singularities of B. Finally Section U] contains the 
computation of the equations of the singular quartic surfaces. We give details 
for two cases, being the others analogous. 

Notation 

By cMTOe we mean algebraic curve over the complex numbers. An (mi,m2, . . .)- 
point, or point of type (mi, m2, . . .), is a singular point of multiplicity mi, which 
resolves to a point of multiplicity m2 after one blow-up, etc. A tangent direction 
of a plane curve singularity of multiplicity n is the direction of a line which cuts 
the singularity with intersection number greater than n. 

The rest of the notation is standard in Algebraic Geometry. 
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2 Computation of plane curves 

2.1 Function LinSys 

Let L be a linear system of plane curves. The Magma function LinearSystem 
computes subsystems of L, but only of curves with ordinary singularities. To 
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overcome this restriction we define the function LinSys{L,p,m,t) (see the Ap- 
pendix), which calculates the subsystem containing the sections which have a 
point of type m[i] = {m[i]i, . . . ,m[i]j) at p[i] with tangent directions given by 
the sequence t[i] of vectors, i — 1, • . • , i^p. Basically this function computes the 
necessary blow-ups and uses the Magma function Linear System after each one. 
Then it blow-downs to return to P^. 

In this function the blow-up at a point with coordinates (x, y) = (a, 6) is 
given by evaluating the elements of L at (x, (x — a)y -|- &), except if the tangent 
direction of the singularity is the one of the vector (0, 1) when it is given by 
evaluating at {{y — h)x + a, y). So, by default, the exceptional divisor correspond- 
ing to the blow-up at (a, b) is the line of equation x = a, except if the tangent 
direction is (0, 1), when it is the line of equation y = b. 

For example, a linear system of cubic curves with an assigned cusp can be 
obtained as follows: 

A<x,y>:=Aff ineSpace(Rationals() ,2) ; 



Let J be a linear system of plane curves. We want to compute the set of tuples 

{pi, . . . ,Pn) such that there is at least an element of J with given singularities 
at the points pi, . . . ,Pn- Thus we have to impose conditions, on generic points, 
to the elements of J. These conditions are minors of a matrix. 

For example, suppose we want to find the elements of J which have a tacnode 
(i.e. a (2, 2)-point). Let (a, b), {u, v) be the coordinates of points pi,P2 such that 
P2 is infinitely near to pi. We suppose that the chart of the blow-up at pi is 
such that the corresponding exceptional curve has equation x — a. Notice that 
then a = u. Denote the equations of the sections of J by ji, . . . ,jm and let 

be the equation obtained from ji by blowing-up at p\, i.e. by evaluating 
ji at (x, (x — a)y -\- b). There are unique polynomials hi,gi, ki such that = 
hi{x — a)^ + gi{x — a) + ki and ki are not divisible by x — a. If the curves 
of J have a double point at pi, then ki = gi = 0, Vi. To obtain double points 
at Pi,P2 we impose conditions, vanishing of derivatives, to the equations ji and 
hi, i = 1, . . . , m. Consider the matrix 



LinSys(LinearSystem(A,3) , A! [0, 0] , [2, 1 , 1] , [ [1 , 1] , [0, 1] ] ) ; 





Resolution of the cusp 



2.2 Procedure CndMt 



Mt := 



.7i(a,&) 
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If there is a tacnode at pi , the maximal minors of this matrix vanish. 

To define such a matrix of conditions in a more general case, involving curves 
with more complicated singularities, we give in the Appendix the Magma pro- 
cedure CndMt{R, J, P, Af, T, ~E, ~Mt). 

Here i? is a polynomial ring, J is a linear system of plane curves, P is 
a sequence of points, an element of M is a sequence of multiplicities and an 
element of T is a sequence of vectors, corresponding to tangent directions of 
infinitely near singularities at an element of P. 

As output we have the matrix Mt and the sequence of polynomials E. In the 
function given in Section r2.3l below. we obtain infinitely near points by imposing 
E = Q. 

2.3 Function ParSch 

The variety S of tuples {pi, . . . ,pn) such that there is an element of J (possibly 
non- reduced) with given singularities at pi , . . . , p„ is defined by: 

• the vanishing of the maximal minors of Mt] 

■ E = (infinitely near points). 

If we are able to compute a point in S, then we use the function LinSys, 
described in Section 12.11 to compute the sections of J which have the given 
singularities at pi, . . . 

In the Appendix we define the Magma function ParSch{J, P, M, T, Eq, Ne, d) 
which computes such a scheme S. Here J is a sequence of linear systems and 
the number d means that the last d linear systems of J do not have certain 
singularities at the points in P. In this case some of the defining equations of S 
are of the type 

{1 + nai) ■ ■ ■ {1 + naa) =0, 

where each is a maximal minor of a certain matrix. This is useful to obtain 
reduced curves. 

The input iVe is used to obtain no repetition on the points. For example, 
suppose we have points pi, . . . ,pi (possibly infinitely near) and we want pi ^ p^ 
and p2 7^ P4- Then we use Ne ~ [[1, 3], [2, 4]]. Analogously Eq is used to obtain 
repetition. 

As in the function LinSys, the exceptional divisor corresponding to the 
blow-up at (a, b) is the line of equation x — a, except if the tangent direction is 
(0, 1), when it is the line of equation y = b. 

3 Quartic surfaces 
3.1 Quartic double planes 

Let A" C P'^ be a normal quartic surface with rational double points as its only 
singularities. Choosing coordinates {x, y, z, w) such that X has a double point 
p at [0:0:0:1], the surface is given by an homogeneous equation 



a2W^ + 63W + C4 = 0, 
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where 02, &3, C4 are polynomials in (a:, y, z). The projection of X from p to is 
a 2 : 1 map, ramified over the sextic curve B with equation 63 — 4a2C4 = 0. The 
image of p in is the conic C with equation a2 — 0, which is tangent to the 
curve B. 

Let X' be the double cover of P^ with branch locus B, i.e. the surface given 
by = 63 — 4a2C4 in the weighted projective space P[3, 1, 1, 1], with coordinates 
{t, X, y, z). The pullback of each irreducible component of C in X' is the union 
of two (— 2)-curves. One can obtain X contracting one of these (— 2)-curves, for 
each component. 

Thus, given such curves B and C, to obtain the equation of X it suffices to 
factor the defining polynomial of B has 63 — 4a2C4. This is not a difficult task. 
A Magma function which computes this factorization and gives the equation of 
X can be found at http://home.utad.pt/~crito/ . 

3.2 Singularities of double coverings 

Each singularity g of X is an A-D-E surface singularity. The corresponding 
singularity of the branch locus B is an A-D-E singularity of the same type, 
except if q is the projection point p f |BHPV[ III. 7.]). In this case the branch 
curve B is tangent to the conic with equation 02 = 0. A generic projection gives 
rise to a singularity in B as in Table [1] (column [q = p]) and conversely one can 
verify that the canonical resolution of B gives graphs as indicated in the table. 
Here the data m (multiplicities) and t (tangent directions) are input data for 
the Magma functions LinSys and ParSch. 

4 Equations 

In this section we compute equations of singular quartic surfaces for the cases 
referred in the Introduction. This is done by first computing the equation of 
the sextic branch curve B. In each case the first singularity corresponds to the 
projection point. For example, in a case E„EiDj we compute a sextic curve with 
singular points as in Table [TJ row E„, column [p ~ q] and rows E^, Dj, column 

4.1 D5E7EV 

In this case B = Bi + B2 is the union of two cubics containing points pi , . . . , p4 
such that: 

Bi , i?2 are tangent to the line piP2 at pi ; 
Bi has a cusp at p^ tangent to B2', 
B2 has a cusp at p^ tangent to Bi ; 
Pi , p2 are not in the line p^Pi 

(without this last condition we would obtain a non-reduced curve). 

We want to find 9 points Pi, p'l, P2, Ps, p'3, p'3 and p4,p'^,p'l with p'^ infinitely 
near to pi and infinitely near to p[. We fix all points except p^ and p'^. 
Notice that each infinitely near point is defined by the tangent direction of the 
corresponding singularity. 

A<x,y>:=Aff ineSpace(Rationals() ,2) ; 

L : = [LinearSystem(A,3) ,LinearSystem(A,3) ,LinearSystem(A, 1) , \ 
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q=p 


Ai 
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\y...)>^...\^.. 


pK2ri — 1 

n > 2 


m=[2^_^ 

n times 


m=[2^_^ 

.^'^ n—1 times 


A2 


y m=[2,l,l] 
<^ [0,1]] 


y J y J y / 


A9„ 

n > 2 


m = [2^_^, 1, 1] 

n times 

^ *=[[],•••, [0,1]] 


X J-' m=[2_^,l,l] 

n—1 times 

/ * = [[],..., [0,1]] 


VjOvi J-A 

n > 


> 


y m=[3,2__^] 

\ n times 


X m=[3, 2^_^] 

n— 1 times 






y m=[3,l,l] 
i=[[],[0,l]] 




D2n+5 
n > 1 




m=[3,2^_^,l,l] 

■^C ri times 

^ *=[[],-•■, [0,1]] 


rr? = f3 2 2 1 ll 

III — i^o, ^, • y , , J-, J-J 

n—1 times 

^ *=[[],••■, [0,1]] 




V / m= [3, 1,1,1] 
*=[[], [0,1], [0,1]] 


....V^.... m=[2,2,2] 

i.n. =3 + 3 


Et 


y m=[3,2,l] 

^ ^ = [[],[o,i]] 


..^^rj,^.. m=[3,2,l] 

i.n. = 1 + 2 + 3 


Eg 


/ m= [3,2,1,1] 
t=[[],[0,l],[l,0]] 


m= [3,2,1] 
...^V... t=[[],[0,l]] 
^3^ *=[[]'[1'0]] 



Table 1: Singularities of branch curve 

exceptional line with self-intersection — 1 

the image of {a2 = 0} in 

i.n. = 'intersection number' 
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LinearSystem(A, 1)] ; 

P: = [A! [0,0] ,A! [1,0] ,A! [2,1]] ; 

M: = [ 

[[1,1] , [1] , [2,1,1] , [1,1,0]] ,\ 

[[1,1],[1],[1,1,0],[2,1,1]],\ 

[[0,0],[1],[1,0,0],[1,0,0]],\ 

[[1,0] , [0] , [1,0,0] , [1,0,0]]] ; 
T: = [[[l, 0]], [],[[!, 2], [0,1]],[[], [0,1]]]; 
S:=ParSch(L,P,M,T, [] , [[1,3,4,7]] ,2) ; 

With the input [1,3,4,7] we obtain no repetition among the points in positions 
1, 3, 4 and 7 (i.e. _pi, . . . ,_P4), so we get p4 ^ pi,p2,P3. The previous sequence [ 
] is empty because we do not want any points to be equal. 

Now it only remains to find a solution in S. Since S has dimension 1, we 
define a zero-dimensional subscheme Si and compute its points. 

R:=Ambient(S) ; 

SI :=Scheme(S, [R.2-2] ) ; 
PoiiitsOverSplittingField(Sl) ; 

The first four coordinates of the ambient space R oi S correspond to p4 and p'^. 
We choose a solution and compute the curve B with the given singularities at 

Pl,...,P4. 

P: = [A! [0,0] ,A! [1,0] ,A! [2,1] ,A! [22/7,2]] ; 
M: = [[2,2] , [2] , [3,2,1] , [3,2,1]] ; 

T: = [[[1,0]] , [] , [[1,2] , [0,1]] , [[1,28/23] , [0,1]]] ; 
J:=LinSys(LinearSystein(A,6) ,P,M,T) ; 

This linear system J has only one section, which factors as 

(x"3 - 323/63*x"2*y - x"2 + 512/63*x*y"2 + 
92/21*x*y - 254/63*y-3 - 88/21*y-2 - 2/7*y) 
(x''3 - 515/126*x~2*y - x'2 + 2482/441*x*y'2 + 
58/21*x*y - 1129/441*y"3 - 317/147*y"2 + 2/7*y) 

Prom this branch curve we obtain the equation of the quartic surface X with a 
D5 point and two E7 points: 

w"2*y"2 + w*x"3 - 129/28*w*x~2*y - w*x~2*z + 25/7*w*x*y*z 

- 856343/254016*x-4 + 1907707/111132*x-3*y - 348881/18522*x-2*y"2 
+ 302119/27783*x*y-3 - 143383/55566*y'4 + 50963/10584*x-3*z 

- 126379/6174*x'2*y*z + 48976/3087*x*y-2*z - 89935/18522*y"3*z 

- 793/588*x'2*z'2 + 5224/1029*x*y*z"2 - 4433/2058*y"2*z"2 

- 17/147*x*z"3 + 299/2058*y*z"3 + l/49*z~4 

4.2 EgDia 

Here B is an irreducible scxtic ciu'vc with a (2, 2, 2)-point and a (3, 2, 2, 2, 2)- 
point which resolves to a cusp after 4 blow-ups. To obtain a reduced ciu-vc we 
impose that B is not a cubic with multiplicity 2 nor a conic with multiplicity 3. 
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A<x,y>:=Aff ineSpace(Rationals() ,2) ; 

L: = [LinearSystem(A,6) ,LinearSystem(A,3) ,LinearSystem(A,2)] ; 

P:=[A! [0,0] ,A! [0,1]] ; 

M: = [[[2,2,2] , [3,2,2,2,2,1,1]] ,\ 

[[1,1,1] , [2,1,1,1,1,0,0]] ,\ 

[[1,1,0] , [1,1,1,1,1,1,0]]] ; 
T: = [[[1,0], [1,0]], [[1,0], [1,1], [],[],[], [0,1]]]; 
S:=ParSch(L,P,M,T, [] , [] ,2) ; 

Now wc proceed as in Section [4. II The equation of the quartic surface X with 
an Eg point and a D13 point is 

w"2*y"2 + w*x'3 - 9/256*w*x~2*y + 69/64*w*x*y*z + lll/64*w*y*z~2 

+ 141393/262144*x-4 - 231/256*x-3*y - 33/128*x-2*y-2 

+ 489/512*x*y~3 - 25/32*y~4 + 57363/32768*x'3*z + 33/64*x~2*y*z 

- 729/256*x*y~2*z + 99/32*y'3*z + 75/32768*x'2*z'2 

+ 1449/512*x*y*z~2 - 147/32*y~2*z~2 - 21/8192*x*z'3 + 97/32*y*z"3 
+ 33/16384*z-4 

4.3 The remaining equations 

Now we give the remaining equations of the quartic surfaces X with rational 
double points with total Milnor number 19 which have no point of type A„. The 
details can be found at http://home.utad.pt/~crito/ . 

w"2*y"2 + w*x"3 - ll/4*w*x~2*y - w*x'2*z + 4*w*x*y*z - 3*w*y*z~2 
+ 113/64*x~4 - 8*x"3*y + 16*x"2*y"2 - 7/2*x*y'3 - ll/2*y~4 

- 6*x~3*z + 25/2*x~2*y*z - 30*x*y~2*z + 23/2*y~3*z + 89/8*x"2*z"2 

- ll/2*x*y*z~2 + 31/2*y~2*z~2 - 9*x*z~3 - 3/2*y*z~3 + 9/4*z"4 

w"2*y"2 + w*x"3 - 3*w*x~2*y - 10/3*w*x*y*z + 7/3*w*y*z~2 

- 29/12*x~4 + 15/2*x"3*y - 17/3*x~2*y~2 - 4/3*x*y~3 + 4*y'4 
+ 2/3*x~3*z + 14/3*x"2*y*z + 20/3*x*y'2*z - 32/3*y'3*z 

+ 5/18*x~2*z~2 - 28/3*x*y*z~2 + 8*y~2*z~2 + l/9*x*z'3 + l/36*z~4 

EgEgDr ( r~2 - 33/73*r + 9/292 = ) 

w"2*y"2 + w*x"3 - 9/8*w*x~2*y + 3/2*w*x*y*z - 9/2*w*y*z~2 
+ l/256*(-292*r + 207)*x~4 + l/64*(-146*r + 69)*x'3*y 

- 15/32*x-2*y-2 + l/64*(146*r - 33)*x*y-3 + l/128*(146*r - 83) 
*y"4 + l/16*(73*r - 66)*x'3*z + 15/8*x~2*y*z + l/32*(-438*r + 99) 
*x*y"2*z + l/16*(-146*r + 75)*y~3*z + 39/32*x"2*z'2 + 1/16 
*(438*r - 99)*x*y*z"2 + l/16*(438*r - 201)*y'2*z'2 + l/4*(-73*r 

+ 3)*x*z~3 + l/4*(-146*r + 59)*y*z~3 + l/16*(292*r - 21)*z"4 

EeEgDs 

w"2*y"2 + w*x~3 - 27/16*w*x~2*y - 21/8*w*x*y*z - 27/16*w*y*z~2 

- 615/1024*x~4 + 27/32*x~3*y - 39/64*x~2*y~2 + 27/32*x*y~3 
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- 25/64*y~4 + 135/256*x'3*z + 39/32*x"2*y*z - 81/32*x*y"2*z 
+ 51/32*y~3*z + 1299/512*x~2*z~2 + 81/32*x*y*z''2 

- 39/16*y"2*z"2 + 351/256*x*z"3 + 53/32*y*z~3 + 297/1024*z"4 

w"2*y"2 + w*x"3 - l/4*w*x~2*y + 3/4*w*x*y*z + l/4*w*y*z~2 

+ 25/64*x"4 - 5/16*x~3*y - l/16*x"2*y"2 + l/8*x*y"3 - l/32*y~4 

+ ll/32*x~3*z + l/8*x~2*y*z - 5/16*x*y'2*z + 3/32*y~3*z 

+ 3/64*x~2*z~2 + l/4*x*y*z~2 - 3/32*y~2*z~2 + l/32*x*z~3 

+ l/32*y*z"3 + l/64*z~4 



E7D 
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w"2*y"2 + w*x"3 - l/2*w*x"2*y + 3/2*w*x*y*z + l/2*w*y*z"2 

+ 13/16*x~4 - 3/4*x'3*y - l/4*x~2*y~2 + 3/4*x*y~3 - l/4*y'4 

+ 5/8*x~3*z + l/2*x''2*y*z - 2*x*y"2*z + 3/4*y~3*z + 3/16*x''2*z''2 

+ 7/4*x*y*z"2 - 3/4*y'2*z'2 - l/8*x*z"3 + l/4*y*z~3 + l/16*z"4 

w"2*y"2 + w*x"3 - 3/4*w*x~2*y - 3/4*w*x*y*z - 3/4*w*y*z~2 

- 23/64*x"4 + 9/16*x~3*y - 3/16*x"2*y"2 - l/16*x*y"3 + l/32*y"4 

- 21/32*x~3*z + 3/8*x~2*y*z + 3/16*x*y'2*z - 3/32*y~3*z 

+ 15/64*x~2*z"2 - 3/16*x*y*z~2 + 3/32*y'2*z"2 + ll/32*x*z"3 

- l/32*y*z"3 + 9/64*z~4 



Appendix: Magma code 

function LinSys(L,p,m,t) 
A : =Ambient (L) ; 

J:=LinearSystem(L,p, [m[i] [1] :i in [l..#m]]); 

for g in [1. .#p] do 

pO: = [BaseField(A) |p[g] [1] ,p[g] [2]] ; 
b: = [];i: = []; 

for j:=l to #t[g] do 

if #Sections(J) eq then break; end if; 

if t[g][j][l] eq then k: = [2,l] ;else k: = [l,2];end if; 

b : =b cat [pO] ; i : =i cat [k] ; 

Bup: = [Evaluate (Sections (J) [i] ,A.k[2] , (A.k [1] -pO [k [1] ] ) *A .k [2] + \ 
p0[k[2]]) div (A.k[l]-pO[k[l]])-m[g] [j] :i in [1 . . #Sections ( J)] ] ; 
p0[k[2]] :=t[g] [j] [k[2]]/t[g] [j] [k[l]] ; 
J : =LinearSystem(LinearSystem(A,Bup) , A!pO,m [g] [j+1] ) ; 
end for; 

for j :=#b to 1 by -1 do 

if #Sections(J) eq then break;end if; 
il:=i[j][l];i2:=i[j][2]; 

Bdn: = [Evaluate((A.il-b[j] [il])"m[g] [j] *Sections(J) [i] ,A.i2, \ 
(A.i2-b[j] [i2])/(A.il-b[j] [il])):i in [1. .#Sections(J)]] ; 
R: =Universe (Bdn) ; 

h : =hom<R->CoordinateRing (A) I [A . 1 , A . 2] > ; 
J:=LinearSystem(A, [h(Bdn[i]) :i in [l..#Bdn]]); 
end for; 



end for; 



return J; 
end function; 

procedure CndMt(R, J,P,M,T,"E,~Mt) 

function D(F,i, j ,a,b) ; 
P:=Parent(F) ; 

for n in [l..a] do F:=Derivative(F,P.i) ;end for; 
for n in [l..b] do F:=Derivative(F,P. j) ;end for; 
return F; 
end function; 

k:=0;a:=0;E: = [] ;Mt: = [[] :i in [1. .#Sections(J)]] ; 
A:=Ambient(J) ; 

PP: = [[P[i]] :i in [l..#P]] cat ID-.i in [l..#M-#P]]; 

for w in [1. .#PP] do 

p:=PP[w] ;m:=M[w] ;t:=T[w] ;tl: = [] ; 

for i:=l to #t do 

if t[l] eq □ then tl: = [t[j] : j in [i . .#t] ] ; t : = [t [j] : j in [l..i-l]]; 

break; end if; 
end for; 

if #in le #t+#p then continue; end if; 

Bup:=Sections(J) ; 

if #p gt and #m gt #t+l then 

pa: = [BaseField(A) |p[l] [1] ,p[l] [2]] ; 
for j :=1 to #t do 

if t[j][l] eq then q: = [2,l] ;else q: = [l,2];end if; 
Bup : = [Evaluate (Hup [i] , A . q [2] , (A . q [1] -pa [q [1] ] ) *A . q [2] +pa [q [2] ] ) \ 
div (A.q[l]-pa[q[l]])-m[j] :i in [l..#Bup]]; 
pa[q[2]] :=t[j] [q[2]]/t[j] [q[l]] ; 
end for; 

Bup := [Evaluate (Bup [i] ,A.2, (A . l-pa[l] )*A. 2+pa[2] ) \ 
div (A.l-pa[l])-m[#t+l] :i in [l..#Bup]]; 
end if; 

a:=a+k;k:=2*(#m-#t-#p) ; 

h:=hoiii<PolynomialRing(J)->R| [R. (a+k-1) ,R. (a+k)] >; 
H: = [h(Bup)] ; 

if p eq [] then u:=R. (a+1) ;else u:=p[l] [1] ;end if; 
su:=[u,R. (a+2)] ; 

o:=0; 

for j in [1. .#iii-#t-#p-i] do 
l:=H[j]; 

H[j] := [Evaluate (Evaluate (l[i] ,R- (a+k-1) ,R. (a+2*j-l)) ,R. (a+k) , \ 
R.(a+2*j)):i in [l..#l]]; 
if tl ne [] and tl[j+#p] ne [] 
then 

if tl[j+#p] [o+l] eq 

then o:=Abs(o-l) ;u:=Abs(o-l)*R. (a+2*j-l)+o*R. (a+2*j) ; 
su:=su cat [Abs(o-l)*u,o*u] ; 

else su:=su cat [Abs (o-l) *u+o*tl [j+#p] [1] /tl [j+#p] [2] , \ 
o*u+Abs(o-l)*tl[j+#p] [2]/tl[j+#p] [1]] ; 

end if ; 

else su:=su cat [Abs(o-l)*u+o*R. (a+2*j+l) ,o*u+Abs(o-l)*R. (a+2*j+2)] 

end if ; 
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l: = [Evaluate(l[i] ,R. (a+k-o) , \ 

(R. (a+k-l+o)-R. (a+2*j-l+o))*R. (a+k-o)+R. (a+2*j-o)) :i in [1. .#1]] ; 

for c in [1 . .m [#t+#p+j] ] do 

l: = [(l[i] -Evaluated [i] ,R. (a+k-l+o) ,R. (a+2*j-l+o) ) ) \ 
div (R. (a+k-l+o) -R. (a+2*j-l+o)) :i in [l..#l]]; 

end for; 

H[j+1] :=1; 
end for; 

mt : = [&cat [[D(H[j] [i] ,a+2*j-l,a+2*j ,b,c-b) :b in [0..c], \ 
c in [0. .m[#p+#t+j]-l]] :j in [l..#H]]:i in [l..#H[l]]]; 
mte: = [[Evaluate(mt[i] [o] , [R.i:i in [l..a]] cat su cat \ 
[R.i:i in [a+k+1 . .Rank(R)] ] ) : o in [1 . .#mt [1]] ] : i in [l..#mt]]; 
Mt: = [Mt[i] catmte[i]:i in [l..#mte]]; 
E:=E cat [R. (a+i)-su[i] :i in [l..k]]; 
end for; 

Mt:=Matrix(Mt) ; 

end procedure; 

function EqNe(R,P,M,T,Eq,Ne,d) 

P: = [[P[i]]:i in [!..#?]] cat LD:i in [1. .#M [!]-#?]] ; 

U: = [] ;k: = l; 

for j:=l to #M[1] do 

if P[j] eq □ then p: = [R.k,R. (k+1)] ;k:=k+2; 

else p: = [P[j] [1] [1] ,P[j] [1] [2]] ;end if; 

U:=U cat p; 

o:=0; 

for h:=l to #T[j] do 

if T[j][h] eq □ then U:=U cat [R.k.R. (k+1)] ;k:=k+2; 
else 

if T[j] [h] [o+l] eq then o:=Abs(o-l) ;end if; 
U:=U cat [Abs(o-l)*p[l]+o*T[j] [h] [2-o]/T[j] [h] [l+o] , \ 
o*p[2]+Abs(o-l)*T[j] [h] [2-o]/T[j] [h] [l+o]] ; 
end if; 
end for; 
end for; 

eql: = [] ; 

for i:=l to #Eq do 

eql:=eql cat [U [2*Eq[i] [m] -1] -U [2*Eq[i] [m+1] -1] :iii in [1 . .#Eq[i] -1] ] \ 
cat [U[2*Eq[i] [m]]-U[2*Eq[i] [m+1]] :m in [1 . . #Eq[i] -1] ] ; 

end for; 

neq : = [] ; z : =Rank (R) -d ; 
for i:=l to #Ne do 

for n:=l to #Ke[i] do 
for m:=n+l to #Ne[i] do 

neq:=neq cat [(l+R.z*(U[2*Ke [i] [n] -1] -U[2*Ne [i] [m] -1] ) )* \ 
( 1+R . z* (U [2*Ne [i] [n] ] -U [2*Ne [i] [m] ] ) ) ] ; 
z:=z-l; 
end for; 
end for; 
end for; 

S:=Scheme(Aff ineSpace(R) ,eql cat neq); 



return S; 
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end function; 

f unct ion ParSch (LL , P , MM , T , Eq , Ne , d) 

if Parent (MM) eq Parent ([ [1]] ) then LL : = [LL] ; MM : = [MM] ; end if; 
TO:=T; 

for j :=1 to #T do 

for i:=l to #T[j] do 

if T[j][i] eq [] then TO[j] : = [T[j] [h] :h in [1. .i-1]] ;break;end if; 

end for; 
end for; 

n:=&+[2*(#MM[l] [i]-#TO[i]):i in [1. .#MM[l]]]-2*#P+d; 
for i:=l to #Ne do 

n : =n+#Ne [i] * (#Ne [i] -1) /2 ; 
end for; 

R :=PolynomialRing(BaseField (Ambient (LL[1] ) ) , Integers () !n) ; 

minors : = [[]: i in [l..#LL]]; 
for z in [1. .#LL] do 
L:=LL[z] ;M:=MM[z] ; 

N: = [[M[j] [i] :i In [1. .#TO[j]+l]] : j in [l..#P]]; 
if #P ne then 

J:=LinSys(L,P,N, [TO[i] :i in [1 . .#P]] ) ;else J:=L; 

end if ; 

CndMt(R, J.P.M.T.'ErMt) ; 
minors [z] :=Minors(Mt ,#Sections (J) ) ; 
end for; 

df: = []; 

if d gt then 
z:=Rank(R) ; 

for i:=#LL-d+l to #LL do 

df:=df cat [&* [l+R.z*minors [i] [j] : j in [1 . .#minors [i] ]] ] ; 

z:=z-l; 
end for; 
end if; 

Sl:=EqNe(R,P,MM,T,Eq,Ne,d); 

S2 :=Scheme (Ambient (SI) ,&cat [minors [i] :i in [l..#LL-d]] cat E cat df ) ; 
S:=Intersection(Sl,S2) ; 

return S; 

end function; 
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